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We study the signal for the detection of quasi-stable supersymmetric particle produced in inter- 
actions of cosmogenic neutrinos. We consider energy loss of high energy staus due to photonuclear 
and weak interactions. We show that there are optimal nadir angles for which the stau signal is a 
factor of several hundred larger than muons. We discuss how one could potentially eliminate muon 
background by considering the energy loss of muons in the detector. We also show results for the 
showers produced by weak interactions of staus that reach the detector. 



I. INTRODUCTION 

Ultrahigh energy cosmic neutrinos couldpotentially 
probe physics beyond the Standard Model Interac- 
tions of UHE neutrinos {Ei, > IQ-'^^eV) with nucleons 
probe center of mass energies above 14 TeV. Some frac- 
tion of these neutrinos may produce supersymmetric par- 
ticles or some other exotic particles. These processes are 
suppressed relative to standard model processes, how- 
ever, in some models interesting signals may arise from 
supersymmetric particles with long lifetimes. In most 
SUSY scenarios, particles produced in high energy col- 
lisions decay immediately into the lightest one and are 
thus hard to detect. However in some low scale super- 
symmetric models in which gravitino is the lightest su- 
persymmetric particle (LSP) and R-parity is conserved, 
the next-to- lightest particle (NLSP) is the charged super- 
partner of the right-handed tau, the stau 2\ . Due to its 
weak coupling to the gravitino, the stau is a long-lived 
particle in these models. For the supersymmetric break- 
ing scale, \/F > 5 X 10^ GeV, the long-lived stau could 
travel distances of the order of 10'' km before decaying 
into the gravitino. 

The distance that staus travel before decaying depends 
on the gravitino mass (or equivalently on the supersym- 
metry breaking scale) and the stau mass. Limits on 
the stau mass of about 100 GeV come from its non- 
observation in accelerator experiments [1, 0, H, @ . Re- 
cently it was proposed that staus produced in high en- 
ergy neutrino interactions, where neutrinos originate in 
astrophysical sources, might be detectable in neutrino 
telescopes 0, 

The cross section for the production of staus in 
neutrino-nucleon scattering [71 is several orders of mag- 
nitude smaller than the neutrino charged-current or 
neutral-current cross section Q. However, once pro- 
duced, the long-lived staus have the potential to travel 
through the earth without decaying and thus open up 
a possibility to be detected in neutrino telescopes. The 
long range of staus could potentially compensate for the 
suppression in the production cross section by increasing 
the effective detector volume and therefore enhancing the 
signal. 



The detection of staus depends on the stau lifetime and 
range, so it is important to determine the energy loss as 
it traverses the earth. The details of the range depend in 
part on the supersymmetry breaking and how the quasi- 
stable stau particle is comprised of the SUSY partners 
of the right-handed and left-handed taus. The electro- 
magnetic energy loss has been shown to have the largest 
contribution from photonuclear interactions for stau en- 
ergies between 10^-10^^ GeV, resulting in a range of lO'' 
km.w.e. for masses of the order of a few hundred GeV 
[l0| . Weak interactions may come into play as well. The 
stau range has been shown to be sensitive to the mixing 
angle of right-handed and left-handed staus. When the 
mixing is maximal, weak interactions act to suppress the 
range at energies above ~ 10^ GeV pjj, however, their 
weak interactions have the potential to produce signals in 
neutrino detectors such as the Antarctic Impulse Tran- 
sient Array (ANITA) [H and the Antarctic Ross Iceshelf 
Antenna Neutrino Array (ARIANNA) [l3| . 

The high energies required for stau production lead 
us to focus on the production of staus in interactions of 
cosmogenic neutrinos as they traverse the Earth and/or 
in the detector. These neutrinos originate from cosmic 
ray protons interacting with the cosmic microwave back- 
ground, 

P7(3K) ^ Ntt 

followed by charged pion, muon and neutron decays. This 
flux is guaranteed as cosmic ray fluxes are measured as 
well as the 3K microwave background. We use a con- 
servative cosmogenic neutrino flux evaluated by Engel, 
Seckel and Stanev (ESS) in Ref. [11] ■ They evaluate the 
neutrino flux associated with the measured cosmic ray 
flux by tracing back cosmic ray propagation through the 
background radiation. Depending on the cosmological 
evolution assumed, the overall normalization of this flux 
has an uncertainty of about a factor of four. In addi- 
tion, neutrinos could be produced at the sources of the 
high energy cosmic rays and those are not included in the 
evaluation of ESS neutrino flux. Thus, the ESS neutrino 
flux is a conservative estimate of the cosmogenic flux. 

The cosmogenic neutrino flux, when neutrino flavor os- 
cillations are not included, peaks at high energies, around 
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10® GeV, and thus it is in the energy range where ANITA 
and ARIANNA have very good sensitivity [H, ^ . The 
neutrino flavor ratio for cosmogenic neutrinos deviates 
from the common 2:1 ratio, due to the neutron decay 
contribution to electron neutrinos [l3] ■ This implies that 
one needs to consider three flavor oscillations when con- 
sidering the cosomogenic neutrino flux that arrives at the 
Earth [1^]. 

We consider cosmogenic neutrinos, their propagation, 
stau production, and subsequent energy loss as it tra- 
verses the earth, for a region of parameter space where 
the staus do not decay over the distances required. We 
compare the resulting stau flux when there is no mix- 
ing between right-handed and left-handed stau and when 
there is maximal mixing. We consider muon-like signals 
(charged tracks) produced by staus and its associated 
background. We discuss the potential for eliminating the 
background by measuring the energy loss, which requires 
large volume detectors. Finally we discuss the showers 
produced in the ice due to stau interactions and its back- 
ground from neutrino-induced showers. 

II. CHARACTERISTIC DISTANCE SCALES 
AND STAU PRODUCTION 




FIG. 1: The parameter space of stau mass m and decay pa- 
rameter F^^^ probed with stau decay lengths corresponding 
to a column depth of 10^ cm w.e. The area to the the right of 
each curve is the parameter space probed by demanding that 
stau does not decay as it traverses this column depth. 



The potential for staus to have a long lifetime is an 
essential feature of this analysis. The lifetime of the stau 
depends on the supersymmetry breaking scale {F^^^) and 
stau mass (m?). It is given by 

CT={— ^ ( 10 km. 1 

V 1014 GeV^ J \ nir J ^ ' 

In Fig. m we show the parameter space (rrir, \fF) 
and the region for which staus with certain energy 
(10®, lO^", 10^^ GeV) have long enough hfetimes, 7cr = 
10^ cm w.e., to travel through the Earth without decay- 
ing for a nadir angle of 88°. Different nadir angles, or 
column depths, correspond to different range of these 
parameters. Clearly there is a large parameter space 
in which staus do not decay, but they may interact via 
charged- current interactions, if there is a maximal mix- 
ing between the right-handed and left-handed stau. Here 
we consider the case when stau mass is 150 GeV and the 
supersymmetric breaking scale F^/^ =10^ GeV. 

Staus arc produced in neutrino-nucleon interactions in 
the Earth with the cross sections 0] about three orders of 
magnitude smaller than standard model neutrino weak 
interactions with nucleons [9] . The staus are produced by 
the decays of heavy quarks and sleptons produced in the 
primary neutrino-nucleon interactions. The cross section 
for the stau production is a t-channel production of a 
slepton and a squark via gaugino exchange. This cross 
section depends on the squark mass, gaugino masses and 
the mass of the left-handed slepton. We take the cross 
section for stau production obtained with rriq = 300 GeV, 
= 250 GeV and mj = 250 GeV f^. Stau cross 
sections obtained using different parameters and different 



SUSY models are shown in Ref.[8]. The cross section 
converted into an interaction length, 

C = {Nacj)-^ , (2) 

is shown for neutrino production of staus by the upper 
solid line in Fig. [2l 

Cosmic neutrino fluxes that reach the Earth get at- 
tenuated as neutrinos traverse the Earth toward the de- 
tector, due to weak interactions, primarily the charged- 
current (CC) interactions. Similarly, once created in 
neutrino-nucleon interactions, the stau flux also gets at- 
tenuated when staus interact weakly. The size of these 
effects can be seen in Fig. [21 where we show neutrino 
interaction length Q and the stau interaction length [ll| 
due to their charged-current interactions with isoscalar 
nucleons, as a function of incident energy. This figure 
shows the stau charged current cross section with max- 
imal weak interactions. We also show for reference the 
column depth for nadir angles of 60° and 80°. 

III. LEPTON AND STAU FLUXES 

A. Cosmogenic Flux 

The cosmogenic fluxes that we use as reference input 
fluxes are the standard evolution and strong evolution 
fluxes of Engel, Seckel and Stanev [13]. The fluxes are 
of neutrinos and antineutrinos, with an electron neutrino 
to muon neutrino flavor ratio that depends on energy. 
Part of this energy dependence comes from the eventual 
neutron decay in fry mr^ production. 
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FIG. 2: Interaction lengths for neutrinos and staus. The solid 
curves, from top to bottom, are the interaction lengths for 
neutrinos to produce staus, for stau charged current (CC) in- 
teractions with maximal weak interactions and for neutrino 
CC interactions. The dotted and dashed lines show the col- 
umn depths for nadir angles of 60° and 80° . 



FIG. 3: Cosmogenic neutrino flux with standard evolution 
and with strong evolution p^ . The upper (lower) solid lines 
represent the total neutrino plus antineutrino flux for strong 
(standard) evolution, dashed lines are the Vfi + Vfj, flux without 
including oscillation effects, and the dot-dashed lines repre- 
sent ffj, + flux with oscillation effects included. 



Neutrino oscillations due to non-zero neutrino masses 
modify the neutrino flavor ratio over the cosniological 
distances they travel. Consequently, we use for the cos- 
mogenic fluxes at the Earth (l5| . 
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where F° are the ESS cosmogenic fluxes The cosmo- 
genic tau neutrino flux F^^ — because of the threshold 
for tau production: the dominant process of resonant A 
production ultimately results in only j/g and v^j_ fluxes at 
the point of production. The expressions in Eqs. (3) and 
(4) are written with the assumption of maximal mixing 
between and Vr- For the numerical work below, we 
have taken 9x2 — 33°. 

In Fig. [3l the total neutrino plus antineutrino cosmo- 
genic fluxes are shown with the upper (strong evolution) 
and lower (standard evolution) solid lines. Without neu- 
trino oscillations, the corresponding -f fluxes are 
shown by the dashed lines. Including neutrino oscilla- 



tions, the V,, 



cosmogenic strong and standard evolu- 



tion fluxes are shown by the dot-dashed lines. Neutrino 
oscillations enhance the incident muon (anti-) neutrino 
flux at low energies due to v^. — > i?^ oscillations. 

Finally, as noted above, the cosmogenic neutrino fluxes 
shown in Fig. [3] are considered "conservative." They cor- 
respond to neutrinos produced by cosmic rays interacting 



with the 3K radiation. Neutrinos could also be produced 
in some cosmic accelerators, such as active galactic nuclei 
and gamma ray bursts, but the fluxes have large theoret- 
ical uncertainties We focus on the cosmogenic neu- 
trinos from cosmic ray interactions with the microwave 
background because they are a "guaranteed" source of 
ultrahigh energy neutrinos. 



Muon Flux 



Muons are produced by u^N charged-current inter- 
actions. Electromagnetic energy loss of muons passing 
through matter together with a survival probability that 
depends on lifetime and energy determines the muon flux. 
Electromagnetic energy loss is described by the average 
change in energy per unit column depth X = pz (in terms 
of density p and distance z) via 



dE 
dX 



-(a + P^E) 



(5) 



The radiative energy loss due to bremsstrahlung, pair 
production and photonuclear scattering, characterized by 
/3^, increases with energy from about ~ 4 x 10~^ cm^/g 
for Ef, - 10^ GeV to about ~ 5 - 6 x 10"^ cmVg for 
Ef, - 10^ GeV 113]. Here we take ~ 4 x IQ-^ cmVg, 
which will correspond to the maximum background for 
the stau charged-track signal due to muons. With a ~ 
2x 10~^ GeVcm^/g for the energies considered here, 
dE/dX ~ —Pfj,E. Eq. ^ combined with effects of the 
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finite muon lifetime on the survival probability Pgi 



dPs. 



Ps, 



dX 



CTpE/m 



leads to 



exp 



-/3^(L-X) 



(6) 



(7) 



The muon flux produced by incident muon neutrinos, 
for column depth L is given by 



m,, I 1 e^'^f'^'^"-^) 



X exp 



cTi_,(3^p \ E^ E^ 



(8) 



where we have made some simplifying approximations. 
We have approximated the neutrino- nucleon charged cur- 
rent differential cross section by 



dE' 



(J,^t,{E,)5{E'-0.^E,) 



where E'^ is the energy of the produced muon. This is 
the energy of the muon before it loses energy via electro- 
magnetic interactions. The flux Fy represents the v^ + v^ 
flux, and is the sum of muon and antimuon fluxes. 

We have also approximated the attenuation of the neu- 
trino flux in transit through a column depth X by the 
shadow factor 



S = expi~a!^''ciE^)NAX 



(9) 



The muon fluxes, as a function of muon energy, are 
shown in Figs. S] and [5] with dotted lines for two nadir 
angles, 80° and 88°. Also shown on these figures are the 
stau fluxes, to which we now turn. 



C. Stau Fluxes 

The stau flux depends on a number of inputs. We use 
the squark and slepton masses indicated in Sec. II yield- 
ing the 1/ — > r interaction length shown in Fig. 2. A 
further choice of parameter is the mixing angle between 
weak isospin zero and weak isospin 1/2 scalars, the right- 
handed and left-handed staus, that yields the mass eigen- 
state of the NLSP stau. We consider two limiting cases: 
staus with no weak interactions, and staus with maximal 
weak interactions plj . 

We start with the case of no weak interactions. In this 
case, the evaluation of the stau flux is similar to that of 
the muon flux. In Ref. lOl, we showed that for staus, a 



reasonable parametrization of the energy loss parameter 
f3r is 

Pr bo + bi HE/Eo) , (10) 

6o = 5 X 10~^ cmVg , 

bi = 2.8 X 10^1" cmVg , 
Eo = 10^" GeV . 

With the expression in Eq. pil)l . the stau survival proba- 
bility and the relation between the initial stau energy and 
the final stau energy as a function of distance is modified 
from the constant (3 case, namely [l^ 



E'~{Er) = Eoexp 

Er 



In- 



En- 



-bi{L-X)-. 



(11) 



P,^„{Er,E^) = exp|^^^ 



— (l + ln(£;^/i?o) 



El 



{l + \n{El/Eo) 



11 1 



cr6op \ Er El 



(12) 



X exp 

This leads to a stau flux of 

Fr{Er,L) ~ 2^ dX j di;,e-^"^^(^-)^^F,(£;„Q) 

X NA<J,^r{E,)5{Er-^-E,^^-^) 

X P,^„{Er,El) . (13) 

The prefactor of 2 accounts for the fact that the staus 
appear in pairs, one from each of the chain of decays of 
the initial squark and slepton. All of the neutrino plus 
antineutrino flavors are included in F,^ {E„ , 0) . 

When weak interactions are included, there is an ef- 
fect due to the attenuation of the staus themselves. The 
survival probability is modified and is given by 



dX 



Acff 
(cT pE / rrLr) 



(tN) 



(14) 
(15) 



The solution of the combined equations Eqs. and 
([5]) is done numerically. 

The resulting stau fluxes, for fixed nadir angles of 80°, 
85° and 88° are shown as a function of stau energy in Fig. 
|4] for the ESS standard evolution flux, and in Fig. \5\ for 
the strong evolution. In these figures, the solid lines are 
with minimal weak interactions, while the dashed lines 
are for maximal weak interactions. 
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FIG. 4: For fixed nadir angles of 80° and 88°, the stau flux 
assuming no weak interactions of the staus (solid lines) and 
assuming maximal weak interactions (dashed lines) produced 
by the ESS neutrino flux with standard evolution, evaluated 
using Eq. (|13p and an input flux with including oscillations. 
The dotted line shows the neutrino induced muon flux at 80° . 
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FIG. 5: For fixed nadir angles of 80° and 88°, the stau flux 
assuming no weak interactions of the staus (solid lines) and 
maximal weak interactions (dashed lines) produced by the 
ESS neutrino flux with strong evolution, evaluated using Eq. 
(|13p and an input neutrino flux with oscillations. The dotted 
line shows the neutrino induced muon flux at 80°. 



IV. STAU SIGNALS 
A. Stau Charged Tracks 

The fluxes of staus and muons depend on the initial 
neutrino fluxes, the production cross sections, and the 
energy loss parameter f}. As noted in Sec. Ill, for the 
stau signal, the initial neutrino flux includes all neutrino 
flavors while for muons, only the oscillated + flux 
contributes. The neutrino energy required to produce a 
stau or a muon of comparable initial energy along the tra- 
jectory through the earth is about 5 times higher for the 
stau because the mean energy of the resulting stau is ap- 
proximately 1/6 of the incident neutrino energy, in con- 
trast to the muon case vifhere E'^ ~ 0.8E^. This means 
we are probing higher energy neutrinos for staus than for 
muons to produce a quasi-stable particle of a given en- 
ergy. On the other hand, the muons lose more energy in 
transit from the production point to the detection point. 

The production cross section for staus is approximately 
three orders of magnitude smaller than the muon pro- 
duction cross section. Once the staus are produced, they 
lose very little energy as they traverse the earth while the 
muon energy loss is of the order 10^ — 10'^ times greater. 
The stau range can be as high as 10* km.w.e. for vanish- 
ing charged-current interactions or suppressed to about 
10^ km.w.e. for maximal charged-current interactions, 
both higher than the muon range. The neutrino atten- 
uation will also have a large effect on the signals. This 
attenuation acts to deplete more muons than staus since 
the muons must be created very near the detector to be 



seen, whereas the staus can be produced farther away. 

These competing effects account for the large 
stau/muon ratio for specific angles and particle energies. 
In Fig. \6\ we show the ratio of the stau flux to the muon 
flux for the angles 80°, 85°, and 88°, as a function of 
energy, for standard ESS evolution. The solid lines show 
the ratio in the scenario where stau weak interactions 
do not occur. The effect of including maximal charged- 
current interactions of the staus can be seen by the sup- 
pression in the ratio, shown with the dashed lines. We 
note that the suppression in the ratio begins to take effect 
between 10® and 10^ GeV, where the stau range begins 
to be affected by stau charged-current interactions [U\ . 
In Fig. [HI the largest ratio occurs for the nadir angle 
of 80° which represents the largest path length through 
the earth of the trajectories shown. As the nadir angle 
is decreased, there is an enhancement in the ratio until 
it reaches a maximum value and then it drops off when 
the neutrino attenuation begins to be the dominant effect 
for the signal. At lower energies, E:^ ~ 2 x 10® GeV, the 
maximum ratio is 125 (at 70°) when stau maximum weak 
interactions are included and it is 280 when there are no 
stau weak interactions. Furthermore, our signal to back- 
ground ratios increase by about 50% (for energy around 
10^ GeV) if we incorporate energy dependent muon en- 
ergy loss, by increasing value to ~ 6 x 10~® cm^/g. 

Fig. [7| shows the same effects for strong ESS evolution 
where the ratio is somewhat reduced. We show in Figs. 
[5] and [5] the ratio for stau/muon fluxes as a function of 
angle for the two fixed final particle energies 10^ and 10^ 
GeV, for ESS standard and strong evolutions. Again, we 
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FIG. 6: For fixed nadir angles of 80°, 85° and 88°, tfie ratio of 
the stau flux assuming no weak interactions of the staus (solid 
lines) and maximal weak interactions (dashed lines) to the 
muon flux produced by the ESS neutrino flux with standard 
evolution, including oscillations. 



FIG. 8: For fixed energies of lO'^ and 10* GeV, the ratio of 
the stau flux assuming no weak interactions of the staus (solid 
lines) and maximal weak interactions (dashed lines) to the 
muon flux produced by the ESS neutrino flux with standard 
evolution, including oscillations. 
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FIG. 7: For fixed fixed nadir angles of 80°, 85° and 88°, 
the ratio of the stau flux assuming no weak interactions of 
the staus (solid lines) and maximal weak interactions (dashed 
lines) to the muon flux produced by the ESS neutrino flux 
with strong evolution, including neutrino oscillations. 



FIG. 9: For fixed energies of lO'' and 10* GeV, the ratio 
of the stau fiux assuming no weak interactions of the staus 
(solid lines) and maximal weak interactions (dashed lines) to 
the muon flux produced by the ESS neutrino flux with strong 
evolution, including neutrino oscillations. 



see the same enhancement of stau/muon ratio shifted to 
the lower angles shown. The strong ESS evolution shows 
the same behavior as standard evolution, again with a 
suppressed ratio. 

Figures OH] show a significant enhancement of the stau 
flux relative to the muon flux for various energies and 



angles. Extracting the stau signal, however, is quite dif- 
ficult. At issue is the fact that apart from stau decays 
and charged current interactions, the stau moves through 
matter in a way similar to a lower energy muon: the stau 
signal is muon-like. 

In terms of electromagnetic energy loss, one of the is- 
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FIG. 10: The inelasticity weighted differential cross section 
for muons for three muon energies: 10^ GeV (solid lines), 
10'^ GeV (dashed hnes) and lO'^ GeV (dot-dashed lines) for 
pair production, bremsstrahlung and photonuclear energy loss 
processes. 



FIG. 11: The inelasticity weighted differential cross section 
for staus for two stau energies: 10^ GeV (dashed lines) and 
10^ GeV (dot-dashed lines) for pair production and photonu- 
clear energy loss processes. Bremsstrahlung energy losses are 
negligible for staus. 



sues is that the average of the energy loss per unit dis- 
tance at high energy scales with f3E. The electromagnetic 
energy loss parameter (3 for pair production and photonu- 
clear contributions to the energy loss scales as the inverse 
mass of the charged particle. The bremsstrahlung pro- 
cess, important for muon energy loss, is negligible for stau 
energy loss since it scales with the inverse mass squared. 

An important element of the energy loss is the electron 
positron pair production cross section. This cross section 
is important at low values of the inelasticity parameter v 
where 



the ratio of the change in muon or stau energy to the 
initial energy. In Fig. [10] we show the differential cross 
section vdajdv scaled by Na/A for muons scattering in 
iron jl^, |2l[. For ease of viewing the small v region, we 
plot the differential cross section versus 1/v. Because the 
photonuclear process includes inelastic scattering correc- 
tions [13, HI], the differential cross section depends on 
the incident muon energy. Pair production dominates 
for \/v greater than ~ 10, that is to say for muon energy 
losses of less than ~ 10% of the initial muon energy. 

The corresponding differential cross section for staus is 
shown in Fig. 1111 This figure shows that the dominant 
process for energy loss of staus is pair production only 
for energy losses of less than v ^ Similar results 

are obtained for ice. 

Figs. [To] and [Tl] can be interpreted in terms of the 
average energy loss {dE/dX). For a stau of energy E^, 
the average energy loss per unit distance is of order of 
that of a muon with energy E^ ^ \Q~^Er ~ m^/my • E^. 



If the average energy loss per unit distance is the only 
observable quantity, the large stau to muon ratios are 
unobservable because one must compare the high energy 
stau flux with the lower energy, but significantly larger, 
muon flux. Quantitatively, this is shown in Fig. [T^ 
Here we take the ratio of the stau flux at energy E:f = 
'nir/fn^j^ ■ E^ to the muon flux evaluated at E^. We note 
that the ratio is peaked at the same angle as in Fig. 8, 
but the ratio is much smaller. This is because (dE/dX) 
matches, for example, muons with energies of 10^ GeV to 
staus with energies more than two orders of magnitude 
higher, where the stau flux is low. 

The w-dependence of the cross sections in Figs. [10] 
and [11] may open a window for differentiating muons and 
staus in the future. Fig. [TT] shows a marked difference 
between muons and staus for large v (small 1/v). An 
energy deposit of 10% of the initial stau energy of Er — 
10^ GeV is observed in a stau signal interaction, with 
very little in the way of comparable energy losses nearby 
since the large v cross section is so low. With such a 
large energy loss (say, AE = O.lEr), there would be no 
confusing the stau with muon with an energy of E^ = 
'ni^i/TiT'T'Er since E"^ <C AE. On the other hand, a similar 
energy deposit could come for a muon with the same high 
energy, but the profile of the track would be significantly 
different, with more large energy losses along the muon 
track where the multiple pair production energy losses 
contribute. Essential here is the large ratio of stau to 
muon masses. The effect is more difficult to observe when 
the masses are closer, e.g., for taus compared to muons 

In Figs. 1131 and [ni we show the average number of in- 
teractions per unit distance of ice rather than the energy 
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FIG. 12: The ratio of stau flux at energy rescaled by m^/m^ 
as a function of the nadir angle for muon energy of 10^ 
GeV and 10* GeV including stau maximal weak interactions 
(dashed lines) and no weak interactions (solid lines). 



loss per unit distance. This is evaluated using 



FIG. 13: The average number of interactions per km of ice 
for muons. The upper line is for v > 0.01, and the lower line 
is for V > 0.1. 



staus, m=150 GeV 
E„=0.1.0.01 E^j^^ 
photonuclear 



M = 



A 



dv 



da 

dv 



(17) 



for density p and column depth X. Here we have taken 
Vmin = f-o/E ~ 0.1,0.01 for muons and staus. The lower 
line on Figs. fT3l and fT4l are for Vmin = 0.1 and the upper 
lines are for Umin = 0.01. A muon which can deposit 
10% of its initial energy of 10^ GeV proceeds to have 
a significantly larger number of subsequent interactions 
in which 10% of its energy is deposited, as compared to 
the stau of the same energy. The scale of the average 
number of interactions is such that one would need an 
observational coverage of a few tens of kilometer distance 
scales to see the distinction between muons and staus. 



B. Stau Showers in Ice 

Staus reaching the detector could interact in principle 
via weak interactions producing showers. Weak interac- 
tions of staus reduce its range [ill], but provide an addi- 
tional opportunity for its detection via interactions in ice 
which can provide a signal for detectors such as ANITA 
[l^ and ARIANNA [ij]. We discuss here showers pro- 
duced by staus and the background due to neutrino in- 
duced showers. 

We have evaluated stau and neutrino showers for dif- 
ferent nadir angles for energies above 10^ GeV. In the 
shower production, we assume maximal mixing for the 
charged-current interactions for staus. The stau showers 
for a given incident angle are determined by Eq. p3p . 



10' 



[GeV] 



FIG. 14: The average number of interactions per km of ice 
for staus with rrir = 150 GeV. The upper line is for v > 0.01, 
and the lower line is for v > 0.1. 



modified to include the probability to produce showers 
in the ice, 



/ ~ dz' 
'-cc 







The stau interaction length due to charged current inter- 
actions is C}j(j. We use the interaction length for maxi- 
mal weak interactions shown in Fig. [21 The pathlength 
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through the ice, zice, for 9 < 88.56° is given by 



z,^^ = Re cose - ^R%cos^e-2REt + t'^, (18) 

where Re represents the radius of the earth and t the 
average ice thickness, taken to be 2 km. For comparison, 
showers due to neutrinos are found from the equation, 

^^.,shr(£^shr, L) ~ F,{E,^,, L)(l - e-^'"/^&c) . (19) 

We compare the resuhs for the showers due to staus 
and neutrinos for three nadir angles: 80°, 85°, and 88° 
in Figs. [Island [TBI AU of the results in this section were 
obtained using the standard evolution cosmogcnic flux. 

In Fig. [151 the showers due to neutrinos are of the 
order 10'^ — 10^ larger than staus at an energy of 10^ 
GeV. The input neutrino flux for the stau signal will be 
at a higher energy than for the direct neutrino production 
of showers. The neutrino flux falls off as a function of 
energy so the ratio of stau to neutrino showers will be 
suppressed in part due to this effect. In addition, the 
stau flux is suppressed due to the small stau production 
cross section. The probability of producing showers in 
the detector is roughly the same for neutrinos and for 
staus with maximal charged-current interactions, so the 
inclusion of the showers is not sufficient to make up for 
the suppression. 

The shape of the stau induced shower flux changes rel- 
ative to the neutrino induced shower flux, so between 
10^ — 10* GeV there is a peak in the ratio of stau to 
neutrino induced showers. This effect is seen in Fig. 1161 
where the peak in the ratio of stau/neutrino showers ap- 
pears, but as indicated, the ratio is small. It is about 
0.009 for 80°. 

Fig. [T71 shows the stau/neutrino ratio for the two fixed 
shower energies, 10^ and 10* GeV, as a function of nadir 
angle. We can see that the ratio is the largest for 10'' GeV 
at about 78°. The integrated flux taking into account the 
contribution due to all incident angles is given in Fig. 1181 
The maximum stau/neutrino ratio occurs for an energy 
of 2.5 X 10'' GeV and corresponds to a stau signal only 
0.33 % of the neutrino signal. 

Whereas muon-like signals from staus can be large 
compared to muon signals from neutrinos because of the 
long lifetime of the stau, we find that the shower signal to 
background ratio is quite small for all angles. This is due 
to combination of the small production cross section for 
stau and the energy effects due to the chain of produc- 
tion and the energy loss of the stau. Stau attenuation 
as it traverses the earth is roughly of the same size as 
neutrino attenuation, as are the probabilities for staus or 
neutrinos to produce showers. 



V. CONCLUSIONS 

We have studied signals of staus produced in interac- 
tions of cosmogenic neutrinos. We have considered two 




Energy (GeV) 



FIG. 15; The shower flux for the incident angles 80° (lower 
curves), 85° (middle curves), and 88° (upper curves) from 
staus (solid lines) and neutrinos (dashed lines). 
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FIG. 16: The ratio of stau to neutrino induced shower fluxes 
for the incident angles 80°, 85°, and 88°, represented by 
curves with the highest, intermediate and lowest peaks, re- 
spectively. 



types of signals, muon-like charged tracks and showers. 
We have focused on low scale supersymmetric models 
that have stau as NLSP, which decays into the lightest 
SUSY particle, the gravitino. For a sufficiently large scale 
of supersymmetry breaking, the stau has a very long life- 
time. 

Our focus has been on cosmogenic neutrino fiuxes and 
their associated stau production in the Earth. Energy 
losses, both through electromagnetic and weak interac- 
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FIG. 17: The ratio of the fluxes of showers due to stau and 
neutrino fluxes for the energies 10^ and 10* GeV as a function 
of incident angle. 
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FIG. 18: The ratio of stau induced to neutrino induced shower 
fliixes, integrated over all incident angles. 



tions, are important in evaluating stau signals. The en- 
ergy loss of staus, however, is relatively small in com- 
parison with muons. Thus, for some nadir angles, the 
stau flux is much larger than the muon flux produced in 
neutrino charged-current interactions. 



The enhancement of the stau flux is larger from an 
input cosmogenic neutrino flux than for the Waxman- 
Bahcall neutrino flux This is because the cosmo- 

genic neutrino flux is peaked at energies of about 10^ 
GeV, while the WB flux is characterized by a steep power 
law with index of two. The large ratio of staus to muons 
from cosmogenic neutrinos is encouraging for experimen- 
tal detection, but in order to see this signal one needs to 
be able to distinguish between staus and muons. Using 
the average energy loss per unit distance is not a good 
way to distinguish staus and muons, since the scaling of 
the energy loss parameter /3 has the effect of making a 
high-energy stau look like a lower energy muon. We have 
proposed a way to distinguish between stau and muon 
tracks by measuring the energy loss of muons via their 
interactions in the ice, and to use this method to reduce 
the background. 

We also considered showers produced by staus inter- 
acting in the ice via charged-current interactions. The 
backgrounds for this signal are showers induced directly 
by neutrinos that reach the detector and interact inside 
the detector via charged-current or neutral-current inter- 
actions. The only way that staus would produce showers 
in the ice is if there is a weak mixing, but this process 
also contributes to reducing stau range. These effects 
combine with the small stau production probability to 
give fluxes of attenuated staus that are several orders 
of magnitude less than attenuated neutrino fluxes. This 
small stau to neutrino ratio translates directly to shower 
rates. 

In addition to weak interactions, another possibility for 
shower production would be stau decays in the detector. 
For the parameter space considered here, with long-lived 
staus, the signal from decays is suppressed relative to 
their weak interactions. 

In conclusion, stau signals at high energies are best 
identified by muon-like tracks. The most important de- 
tection issue is distinguishing between staus and muons, 
which may be possible by looking at the incremental elec- 
tromagnetic energy loss as the charged particle moves 
through the detection volume. With very large volumes, 
there is a potential for detection of staus with future large 
neutrino telescopes. 
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